On the sign structure of doped Mott insulators 
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We demonstrate that the sign structure of the t- J model on a hypercubic lattice is entirely different 
from that of a Fermi gas, by inspecting the high temperature expansion of the partition function 
up to all orders, as well as the multi-hole propagator of the half-filled state and the perturbative 
expansion of the ground state energy. We show that while the fermion signs can be completely 
gauged away by a Marshall sign transformation at half-filling, the bulk of the signs can be also 
gauged away in a doped case, leaving behind a rarified "irreducible" sign structure that can be 
enumerated easily by counting exchanges of holes with themselves and spins on their real space 
paths. Such a sparse sign structure implies a mutual statistics for the quantum states of the doped 
Mott insulator. 

PACS numbers: 71.10.Fd,71.27.+a,02.70.Ss,74.20.Mn 



The progress in the understanding of the physics of 
strongly interacting electron systems has been strongly 
hindered by the infamous fermion minus sign problem 
rendering field theoretical and statistical physics meth- 
ods to be ill behaved for fermions. The t-J model, 
catching the essence of the doped Mott insulators, is 
archetypical. Despite twenty years of concerted effort, 
inspired by its relevance towards the problem of high Tc 
superconductivity, f?| nothing is known rigorously about 
this model, except then for the one dimensional case. In 
fact, the other exception is the Mott-insulating state at 
half-filling, where the Hubbard projection turns the indis- 
tinguishable fermions into distinguishable spins, and the 
remnant signs of the unfrustrated spin problem can be 
gauged away by a Marshall sign transformation. [^I Upon 
doping, however, the fermion signs get active again but it 
is obvious that the sign structure has to be quite different 
from that of a Fermi gas, given that all signs disappear 
at half filling. 

It is instructive to first specify the sign structure in a 
Fermi gas. In a path- integral formalism, 3] the partition 
function of a Fermi gas can be expressed as 
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with each path c composed of a set of closed loops of 
the spatial trajectories of all fermions and Zq [c] > 0. 
The sign structure is then governed by (— l)^"^!"^!, with 
Ncx[c] = N — N\oop[c\ where N = J2w wCw{c) is the total 
number of fermions and iVioop [c] = J2w ^™ ^^"^ closed 
loop number, in which w denotes the number of fermions 
in a loop (also called the winding number Q of the loop) 
and Cw (c) , the number of loops with a given w for a given 
path c. 

Here we report our discovery of a remarkably sparse 
sign structure for the t-J model, which can be rigorously 
identified at arbitrary doping. Basically, we shall prove 



that the partition function for the t-J model is given by 



where -E[c] > [see ([T4|) ] and the sign structure 
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for a given c composed of a set of closed loops for all 
holes and spins (an example is shown in Fig. [T]), where 
[c] denotes the total number of exchanges between the 
holes and down spins and -/V^xl^] the total number of ex- 
changes between holes. In addition to appearing in the 
above partition function, the sign structure Tc will be 
also present in various physical quantities based on ex- 
pansions in terms of quantum-paths in real space: the 
n-hole propagator of the Mott-insulating state as well as 
the zero temperature perturbation theory of the ground 
state energy (both up to all orders). 

Compared to the full fermion signs in ([T|), which is 
an exactly solvable problem for a Fermi gas,[3| the "sign 
problem" for the t-J model then becomes that Tc in is 
too sparse to be treated as a fermion perturbative prob- 
lem. It implies that in the mathematically equivalent 
slave-boson representation, the no double occupancy con- 
straint must play a crucial role to "rarefy" the statistical 
signs of fermionic "spinous" in order to reproduce the 
correct sign structure Tc, which disappears at half- filling. 
On the other hand, in the slave-fermion representation 
besides the statistical signs associated with the fermionic 
"holons" [related to in ([3])], extra signs in Tc will have 
to be generated dynamically, which are previously known 
as the phase strings identified in the one-hole case.[l| 

In particular, we will show that in the two-dimensional 
(2D) case Tc can be precisely captured by a pair of mu- 
tual Chern-Simons gauge fields: the electrical charges 
feel TT flux-tubes attached to the spin "particles" and 
vice versa, in an all-boson formalism which is known as 
the phase string formulation derived before by a different 
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FIG. 1: A typical diagram for a set of closed paths, denoted by 
c in the expansion of the partition function given in (|2]) . Only 
the hole and down spin loops are shown as the up spins are 
not independent due to the no double occupancy constraint. 
Pure spin loops without involving exchanges with the holes do 
not contribute to any signs and are not explicitly shown here. 
The total sign Tc associated with the diagram, defined in ([Sjl , 
is determined by counting the hole-down-spin exchanges and 
the hole permutations. In this particular c, A^^ = 1 + 2 + 1 = 4 
and iV4 = + 2 + 1 = 3 such that Tc = -1. 

method. So the unusual sign structure Tc strongly hints 
a mutual statistics nature of this doped Mott insulator, 
and thus offers critical guidance in the construction of 
correct quantum states of it. 

Let us begin with the t-J model on a bipartite lattice 
of any dimensions Ht-j = Ht + Hj. In the slave-fermion 
representation, the electron annihilation operator can be 
written as 

C^a = {-oYflb,, (4) 

where / denotes the fermionic holon operator and h the 
bosonic spinon operator, which satisfy the no double oc- 
cupancy constraint ///^ + ^L^icr = 1. Then the hop- 
ping and superexchange terms can be expressed, respec- 
tively, as follows 



where 
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Here Poi and Poi denote the nearest neighbor hole-spin 
exchange operators, the nearest neighbor spin-spin 



exchange operator, while the Q term describes the po- 
tential energy between the nearest neighbor antiparallel 
spins. 

Note that the Marshall signQ factor (— cr)' is explic- 
itly introduced in ^ such that the superexchange term 
Hj acquires a total negative sign in front of the spin 
exchange and potential operators. Then one finds the 
matrix element {(j>'\Hj\(j)) < where |0) and denote 
the Ising spin basis bl^^^bl^^^...\Q), which imphes that the 
Hj term will not cause any sign problem. In particular, 
the ground state of Hj can be always written as 

which is true even at doped case so long as there is no 
hopping term. 

Partition function. The nontrivial sign problem only 
arises when holes are doped into the system and allowed 
to hop. It can be traced to the sign difference between the 
hole-spin exchange operators. Pot Sind Poi, in the hop- 
ping term ([5|) in addition to the sign problem associated 
with fermionic holons. By making the high-temperature 
series expansion of the partition function up to all orders 

Zt^j = Trle-/'^*-^}^ Vt^Tr{(i/,_,;n 
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and inserting the complete set 

E = 1 (13) 

<t>{ih} 

between the operators inside the trace (here \(f); {Ih}) is an 
Ising basis with specifying the spin configuration and 
{Ih} denoting the positions of holes), one can evaluate 
term by term of the expansion in ([12]). Because of the 
trace, the initial and final hole and spin configurations 
should be the same such that all contributions to Zt-j 
can be characterized by closed loops of holes and spins 
although each of them can involve multi-holes or -spins 
as shown in Fig. [1] 

Finally one arrives at the compact form given in Q 
based on the above high-temperature expansion, with 

/2A^-^-'-' (/3J/2)» 

2^ -^^ SnM,.+Mu+MQ (14) 

^ ^ n 

in which A//i[c] and M||[c] represent the total steps of 
the hole and down spin "hoppings" along the closed loops 
for a given path c, and Mq [c] the total number of down 
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spins interacting with the nearest-neighbor up spins via 
the potential term Q in ([6]). Obviously Z[c] > in p4)) . 
Thus the nontrivial sign structure of the partition func- 
tion Zt^j in ^ is entirely captured by in ([3]) where 
iV^ [c] denotes the total number of exchanges between the 
holes and down spins [i.e., those actions taken via Poi in 
([5|)] and iV^xfc] the total number of exchanges between 
holes arising from the fermionic statistics of the holon 
operator /. 

The expression ^ for the partition function clearly 
demonstrates that precisely depicts the irreducible sign 
structure at arbitrary doping, temperature, and dimen- 
sions for the t-J model on a bipartite lattice. In the 
following, we shall further illustrate how Tc similarly ap- 
pears in other physical quantities. 

Multi-hole propagator. Define the multi-hole propaga- 
tor 

Gi{js},{ishE) 

= ii'olcl^^.cl^^, ■ ■ ■ G{E) ■ ■ ■ C^^^,Ci,aAi'0} (15) 

where \^o) is the half-filling ground state and 
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One can make the following expansion which converges 
at £^ < Eg (the multi-hole ground state energy): 
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and then insert the complete set p3|) between the ex- 
change operators. Similar to the evaluation of the par- 
tition function, denoting c as a given set of open paths 
connecting the hole configurations {is} and {js}, with 
IV'o) expanded in terms of |0) [ PT]) ]. we find 

G({j.}, {^s}■, E) = ~AY^ ^ ^^^i^' ^] (18) 

in which each set of paths c is weighed by the phase 
strings Tc and an amplitude 



W[c;E]=[^ 
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with A = (-o-J*°~^° . At E < Eg < 0, the expan- 

sion (fT8)) is converged, and W[c; > shows that Tc 
is indeed an "irreparable" (irreducible) sign which is ex- 
pected to play a critical role via constructive and destruc- 
tive quantum phase interferences among different "path" 
c's. Note that the single-hole version of ([T5)) has been 
previously discussed in Ref . 0, 0] . 



Ground state wave function. Define a wave function 
vl/o[7^]^(V'o|c,U,<,^•••|*o) (20) 

with I^Po) as the true ground state and TZ = {ih}', Ws}- 
Then, according to (fT5|) and (fT6|) , [Ti] will be selected 
as i? ^ Eg from below, with (|18p implying 



(21) 



where on the right hand side the path c^'s are all the 
closed loops connected to TZ, each weighed by a positive 

amplitude W[cn] = X^X^'^^W. E] 



-Bo 



E^Eo 



Therefore the sign structure Tc-j^ must be naturally built 
into the ground state wave function. In the following 
we first examine the ground state energy, and then the 
mutual statistics implied for the wave function. 

Ground state energy. Based on the Goldstone's 
theorem,!^ the energy shift of the ground state due to 
the hopping term Ht can be expressed by 



oo 

E 
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1^ 



/connected 



(22) 

where -ffj|\E'o) = fiol^o) and the subscript "connected" 
means that only matrix elements of the operator in ([22]) 
which start from the ground state j^'o) and end with I^Pq) 
without disconnected parts should be included. 

Here j^'o) is generally written in a translational invari- 
ant form with a momentum K: 



l*o(K)) 



R 



where |<I>o; {r/,^ }) is also the ground state of Hj for a 
set of hole distribution {r/^} which minimizes the su- 
perexchange energy fig. Like in the half- filling case, 
one can expand |<I>o; {r/,^ }) in terms of the Ising basis: 
l*o;{rjJ) = E0X0({r/J)l'/';{r;J) with > as 
mentioned before. 

By making the expansion in terms of Hj/Qq and using 
a similar procedure in dealing with the expansion (jl7p . 
one finally gets 

EG-no = noY,e''^'^^~^'^T.x<P'X^ E ^cW[c;no] 

R.R.' (pcj)' c(connected) 

(23) 

in which the path c starts from |<I>o;{r/^ ^ R-}) and 
ends with |$o; {r/i_^— R'}) without including the "dis- 
connected" paths. [6| Indeed the sign factor Tc, weighed 
by VF[c;Oo] > defined in (|19p. determines the ground 
state energy shift upon doping. 

Mutual statistics. Tc in ^ suggests that statistics signs 
of the fermionic holes associated with the slave-fermion 
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representation are actually indistinguishable from the 
"phase strings" generated by the motion of holes. It will 
be thus instructive to treat the whole sign structure on an 
equal footing: take the holon and spinon mathematically 
as all bosons and redefine the hole-spin and spin-spin ex- 
change operators ([5]) and © by 

Po^ = E (^^''''^'^1'^^) (^V't) + H.c. (24) 

(ij) 

-Pol = J2[e-^''^^hlh,) [e-^^^^b]^h^) + H.c.i25) 

Pn - E (^V^t) (e^'^'^'^Va) + H.c. (26) 

(v) 

where the fermionic holon fi is replaced by a bosonic hi 
and the minus sign in front of Pqi is absorbed. The Q 
term pil)) remains unchanged. In 2D case, it is straight- 
forward to verify that if F, , and G,, are chosen as 

F^, = T.im-0M«+O (27) 

G^J = (28) 

where n'l^ and n^'' are the number operators of spinon (cr) 
and holon, respectively, and 9i(l) = Sln(zi — z;) with Zi 
denoting the complex coordinate of site i, then the par- 
tition function ([2]) can be correctly reproduced. Without 
and G^, by contrast, one finds Tc = 1 in ([2]). Namely 
the sign structure is indeed entirely captured by the phase 
factors, e~'^'3 and e~*'^j% in this bosonic formalism. 
Rewriting F,^- = - Af^- + 4 + A^^-, and G,^. = 2^^^., 

and using the constraint X]<t '^/'cr + = 1; o^'^ can 
show that the three link variables, A^^, A'^^, and (/)°-, 

satisfy Er^y- = i'^EieSr ('^n " "u) ' Sr^y" = 
±7rEigEr"('' ^ loop r enclosing an area Er, and 
En = i"" for each plaquette. So they describe tt 
flux tubes bound to spinous, holons, and each plaquette, 
respectively. Since the Hamiltonian is invariant under 
gauge transformations hi /lie"^', Afj Afj + {ipi — ipj) 
and b,a 6*ae*'"^', A'}^ A^^ + (6*, -9,), the bosonic 
holons and spinons carry the gauge charges of Al^ and 
A^^ , respectively. With the effect of Tc described by the 
mutual Chcrn-Simons gauge fields, Af,- and A^-, the t- 
J model in the bosonic formalism explicitly becomes a 
mutual fractional statistics (mutual semions) problem. 
Correspondingly, the electron wave function ipe of the t- 
J model can be also expressed in terms of ipb in this 
bosonic formalism via ■(/'e = A^V'fcjQ i'^ which the large 
gauge transformation /C will transform by 

/C — > TcK. (29) 

under an operation that the hole and spin coordinates 
are continuously permuted via a series of nearest neigh- 



bor exchanges (with the no double occupancy obeyed at 
each step), with the coordinates forming closed loops, 
denoted by c, after the system back to the original con- 
figuration at the end of the operation. By comparison, 
the sign structure in ([1]) is related to the usual antisym- 
metric fermionic wave functions for a Fermi gas. 

In summary, we have demonstrated rigorously that the 
Hubbard projections inherent to the physics of doped 
Mott insulators change the rules of fcrmion statistics fun- 
damentally as compared to the Fermi gas. Pending the 
doping level, the irreducible sign structure that is of rel- 
evance to the physics is much more sparse in the former 
and we have shown that at least in real space expansions 
these irreducible signs are easy to count. In particular, 
in the 2D case, we have established a precise relation in 
which the physical sign structure of the t-J model is ex- 
plicitly determined by the mutual Chern-Simons fields, 
with the wave function satisfying the mutual statistics. 

This does not mean that we have solved the problem - 
the "mutual Chern-Simons" theoryiS] of the phase string 
formulation is still far from being completely understood. 
However, our results open up new alleys for investiga- 
tion. High temperature expansions should be revisited 
to study in detail in what regard the t-J signs differ from 
those of a Fermi gas. It would be quite interesting to find 
out how the "irreducible" hyper nodal surfaces of numer- 
ically determined t~ J model ground states look like. At 
the least, it seems possible to critically test Anderson's 
conjecture[9| that the ground state of the doped Mott- 
insulator has to be orthogonal to that of the Fermi liquid, 
using the elementary fact that wave functions having a 
qualitatively different nodal surface cannot overlap. 
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